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Bone marrow-derived mesenchymal stem cells (MSCs) are 
multipotent, with the ability to differentiate into different cell 
types. Additionally, the immunomodulatory activity of MSCs 
can downregulate inflammatory responses. The use of MSCs 
to repair injured tissues and treat inflammation, including in 
neuroimmune diseases, has been extensively explored. 
Although MSCs have emerged as a promising resource for 
the treatment of neuroimmune diseases, attempts to define 
the molecular properties of MSCs have been limited by the 
heterogeneity of MSC populations. We recently developed a 
new method, the subfractionation culturing method, to iso- 
late homogeneous human clonal MSCs (hcMSCs). The 
hcMSCs were able to differentiate into fat, cartilage, bone, 
neuroglia, and liver cell types. In this study, to better under- 
stand the properties of neurally differentiated MSCs, gene 
expression in highly homogeneous hcMSCs was analyzed. 
Neural differentiation of hcMSCs was induced for 14 days. 
Thereafter, RNA and genomic DNA was isolated and sub- 
jected to microarray analysis and DNA methylation array 
analysis, respectively. We correlated the transcriptome of 
hcMSCs during neural differentiation with the DNA methyl- 
ation status. Here, we describe and discuss the gene ex- 
pression profile of neurally differentiated hcMSCs. These 
findings will expand our understanding of the molecular prop- 
erties of MSCs and contribute to the development of cell 



therapy for neuroimmune diseases. 
[Immune Network 201 4; 14(1 ):54-65] 

INTRODUCTION 

Adult stem cells have lineage-restricted differentiation poten- 
tial. They are generally named according to their tissue of 
origin, for example, as mesenchymal stem cells (MSCs), hem- 
atopoietic stem cells, and endothelial progenitor cells. MSCs 
have been isolated from various mesenchymal tissues, includ- 
ing bone marrow (BM), adipose tissue, and umbilical cord 
blood (1,2). Because of their multifaceted immunomodulatory 
functions, MSCs have emerged as an attractive therapeutic ap- 
proach for the treatment of immune diseases (3). A number 
of saidies have shown that MSCs are effective in treating im- 
mune-related or inflammatory diseases such as graft-ver- 
sus-host disease, colitis, and pancreatitis (4-6). 

Multilineage plasticity is another propeity of MSCs. MSCs 
can differentiate into various cell types found in mesenchymal 
tissues, such as osteocytes, chondrocytes, and adipocytes (7). 
MSCs can differentiate into ectodermal, endodemial, and mes- 
odermal lineages (8,9). MSCs have been used to repair spinal 
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cord injury in vivo and to promote neuronal recovery (10). 
In addition, MSCs can adopt neuronal and glial phenotypes 
in vitro under the appropriate conditions (11,12). Whether 
the neuronal morphology obseived during the neurogenic dif- 
ferentiation of MSCs is an artifact induced by the neurogenic 
medium is controversial (13,14), but several studies have pre- 
sented evidence supporting the neural differentiation poten- 
tial of MSCs (15,16). 

Furthermore, MSCs are considered a potent therapeutic tool 
likely to have practical use in the near future. There have 
been extensive efforts over the past decade to treat various 
incurable diseases using stem cells, including neuroimmune 
disorders such as spinal cord injury, stroke, and multiple scle- 
rosis (17). Because they are multipotent and immunomo- 
dulatory, neural stem cells or neural preairsor cells are a ther- 
apeutic option for the treatment of neuroimmune diseases 
(18,19). However, considering the simplicity of stem cell iso- 
lation, the ease of cell expansion, and the wide range of ap- 
plicability, MSCs offer a good alternative for the treatment of 
neuroimmune diseases, particularly given their neuroglial po- 
tential and immunomodulatory properties. Therefore, a de- 
tailed characterization of neurally differentiated MSCs is need- 
ed to improve cell-based treatments for neuroimmune dis- 
eases. 

The density-gradient centrifugation method is the most 
popular technique for isolating MSCs from BM (20). MSCs ob- 
tained by this method are heterogeneous because they con- 
tain mixed populations of MSCs. Although all of the cells pos- 
sess MSC characteristics, their cell surface marker expression, 
differentiation potential, and cytokine secretion vary, suggest- 
ing that BM contains a variety of MSC populations with differ- 
ent biological capacities (21). Studies of heterogeneous MSCs 
have been informative, but the heterogeneity of the cell pop- 
ulation likely affects the interpretation of data acquired from 
experiments with these cells. In light of this shortcoming, we 
recently developed a novel method, the subfractionation cul- 
turing method (SCM), to isolate and establish homogeneous 
human clonal MSCs (hcMSCs) from small aspirates of human 
BM (21). hcMSC lines estalilished with our protocol express 
well-known MSC markers and differentiate into various cell 
types, including osteocytes, chondrocytes, and adipocytes. 
Additionally, the hcMSC clones express neural or hepatocytic 
phenotypes after neural or hepatogenic differentiation in 
vitro. 

Although several microarray analyses have examined the 
expression profiles of MSCs (22,23), we could not find any 



studies that analyzed gene expression in homogeneous clonal 
MSCs during differentiation. Because a number of MSC pop- 
ulations with different differentiation potentials might be pres- 
ent in the BM (21), we used hcMSCs with high neurogliogenic 
potency to compare gene expression before and after neural 
differentiation. In this study, we verified the neural potency 
of highly homogeneous hcMSCs and used microarray analysis 
to assess gene expression during their differentiation into 
neuroglial cells. Additionally, we performed CpG methylation 
array analysis and compared methylation status with gene ex- 
pression to obtain a more reliable profile of gene expression 
in neurally differentiated hcMSCs. 

MATERIALS AND METHODS 

Isolation of hcMSCs and cell surface phenotyping 
using flow cytometry 

BM aspirates were taken from the iliac crest of a healthy male 
donor after informed consent was provided (approved by the 
INHA University Medical School Institutional Review Board; 
IRB Number 10-51). Isolation of hcMSCs was carried out as 
previously described (21), Several cell surface antigens on the 
established hcMSC line, named KBHD502, were characterized 
by flow cytometry. The antibodies used for the analysis were 
anti-CDl4, anti-CD29, anti-CD31, anti-CD34, anti-CD44, an- 
ti-CD73, anti-CD90, anti-CD105, anti-CD119, anti-CD133, an- 
ti-CDl66, anti-HLA class I, anti-HLA-DR, anti-Stro-1, an- 
ti-c-Met, and anti-c-Kit (BD Biosciences Pharmingen, San 
Diego, CA, USA). The cells were analyzed in a FACSCalibur 
flow cytometer (BD Biosciences). Isotype-matched control 
antibodies were used as controls. 

In vitro immunosuppression assay 

The in vitro immunosuppressive activity of hcMSCs was de- 
termined by [^Hl-thymidine incorporation. Briefly, 2 x lo' pe- 
ripheral blood mononuclear cells (PBMCs) from two different 
healthy donors (1 x lO' cells each) were mixed and cultured 
in a 96-well plate for mixed lymphocyte reactions. hcMSCs 
(4x10^* cells) were co-cultured at a ratio of 1:5 (hcMSCs: 
PBMCs) in these reactions. pH]-thymidine (1 /^Ci/reaction) 
was added for the last 12 ~ I6 h of culture. Radioactivity was 
measured in a beta-counter. 

In vitro neural differentiation and immunofluore- 
scence (IF) staining 

For IF staining, the cells were seeded onto an 8-weIl chamber 
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slide (Nunc, NaperviUe, IL, USA) at a density of 1 x lO"* cells/ 
well. After a 24-h incubation, the growth medium was re- 
moved and replaced with neurogenic dilTerentiation medium 
(neurobasal medium supplemented with B27 supplement 
(Gibco-BRL, Gaithersburg, MD, USA), 1 mM dibutyryl cAMP 
(Sigma-Aldrich, St. Louis, MO, USA), 0.5 mM l-methyl-3-iso- 
butylxanthine (Sigma-Aldrich), 20 ng/ml human epidermal 
growth factor (Sigma-Aldrich), 40 ng/ml basic fibroblast growth 
factor (Sigma-Aldrich), 10 ng/ml fibroblast growth factor 8 
(Peprotech, Rocky HiU, NJ, USA), and 10 ng/ml brain-derived 
neurotrophic factor (R&D Systems, Minneapolis, MN, USA). 
The cells were cultured in serum-free neurogenic differentia- 
tion medium for 2 weeks. At the end of the differentiation 
period, the cells were fixed with 4% paraformaldehyde and 
permeabilized with 0.5% TritonX-lOO/PBS. The cells were la- 
beled with primary antibodies (1 : 200 ~1 : 1,000), including 
rabbit anti-human glial fibrillary acidic protein (GFAP; Sigma- 
Aldrich), mouse anti-neuron-specific class III /? -tubulin (Tujl; 
MiUipore, Billerica, MA, USA), and rabbit anti-microtubule- as- 
sociated protein 2 (MAP2; MiUipore) overnight at 4°C. After 
incubation with the primary antibodies, the cells were in- 
cubated for 1 h with AlexaFluor488- or AlexaFluor594-conju- 
gated secondary antibodies (1 : 300; Molecular Probes, Carls- 
bad, CA, USA). The cells were subsequently stained with 
4' ,6-diamidino-2-phenylindole (DAPI; Molecular Probes) or 
propidium iodide (PI; Molecular Probes) for 1 min. After 
mounting, the samples were analyzed by confocal micro- 
scopy (Zeiss LSM510 Meta Confocal Imaging System; Carl 
Zeiss, Thornwood, NY, USA). In vitro adipogenic, chondro- 
genic, hepatogenic, and osteogenic differentiation are de- 
scribed in the supplemental material (Supplemental Method 
1). 

Isolation of RNA and RT-PCR 

At the end of neurogenic differentiation, total RNA was ex- 
tracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). 
cDNA was synthesized using the Reverse Transcription 
System (Promega, Madison, WI, USA). Semi-quantitative RT- 
PCR was performed using ExTaq DNA polymerase (Takara 
Biotechnology, Shiga, Japan). Amplified PCR products were 
electrophoresed on 1% agarose gels. PCR primer sequences 
are summarized in the supplemental material (Supplemental 
Method 2). 



Microarray and genome-wide CpG methylation 
microarray 

At the end of the neurogenic differentiation period, total RNA 
was isolated using TRIzol reagent. The synthesis of Target 
cRNA probes were synthesized and hybridized using Agilent's 
Low RNA Input Linear Amplification kit (Agilent Technology, 
Palo Alto, CA, USA) according to the manufacturer's instruc- 
tion. The Cy3/5-labeled cRNA amplicon was purified on a 
cRNA Cleanup Module (Agilent Technology). Labeled cRNA 
was quantified using an ND-1000 spectrophotometer 
(NanoDrop Technologies, Wilmington, DE, USA). After chec- 
king the labeling efficiency, the cRNA was fragmented by 
adding blocking agent and fragmentation buffer and incubat- 
ing at 60°C for 30 min. The cRNA fragments were re- 
suspended in hybridization buffer and directly pipetted onto 
Agilent's Whole Human Genome Oligo Microarray (44K). 
The arrays were hybridized at 65°C for 17 h in an Agilent 
Hybridization oven (Agilent Technology). According to the 
manufacturer's protocol, the hybridized microarrays were 
washed. For microarray data acquisition and analysis, the hy- 
bridized arrays were scanned with Agilent's DNA microarray 
scanner and quantified with the Feaaire Extraction software 
(Agilent Technology). Data was normalized, and regulated 
genes were identified with GeneSpring GX 7.3.1 (Agilent 
Technology). The averages of the normalized ratios were cal- 
culated by dividing the average of the normalized signal 
channel intensity by the average of the normalized control 
channel intensity. Functional annotation of genes was carried 
out according to the guidelines of the Gene OntologyTM 
Consortium (http://www.geneontology.org/index.shtml) us- 
ing GeneSpring GX 7.3.1. 

For CpG methylation array analysis, genomic DNA samples 
from control and neurogenically differentiated hcMSCs were 
extracted using a Wizard Genomic DNA Purification kit (Pro- 
mega). Genomic DNA was sonicated to an average fragment 
length of 300 bp using a VCX130 ultrasonic processor (Sonics 
& Materials, Newtown, CT, USA) and purified using a 
QIAquick PCR Purification kit (Qiagen, Valencia, CA, USA). 
To isolate methylated genomic DNA, 2 iig oi His6-tagged 
methyl DNA-binding protein was incubated with 500 ng of 
sheared genomic DNA fragments from control and differ- 
entiated hcMSCs for 8 h at 4°C. The enriched methylated 
DNA was amplified using a Whole Genome Amplification kit 
(Sigma-Aldrich). The amplified products from control and dif- 
ferentiated cells were labeled with Cy3-dUTP and Cy5-dUTP, 
respectively, using the BioPrime Total for Agilent aCGH sys- 
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tem (Invitrogen). The labeled DNA samples were mixed and 
hybridized to a human 244K CpG island microarray (Agilent 
Technology). The slides were scanned with an Agilent scan- 
ner, and the images were quantified with the Feature 
Extraction software (Agilent Technology). Two independent 
experiments were performed. Data were normalized and dif- 
ferentially regulated genes were identified using GeneSpring 
GX 7.3.1 software (Agilent Technology). 



RESULTS 

We isolated and established several hcMSCs from a healthy 
donor's BM by SCM. One clone exhibited excellent multi- 
lineage plasticity. The clonal cells expressed known MSG 
markers, including CD29, CD44, CD73, CD90, CD105, CD133, 
CDI66, HLA-class I, and Stro-l; however, they were negative 
for hematopoietic and endothelial markers, including CD14, 
CD31, CD34, and HLA-DR (Fig. lA). Additionally, the cells 
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Figure 1. Characterization of IncMSCs and neural differentiation. The hcMSC line was established by SCM. (A) The cell surface markers were 
analyzed by flow cytometry. The cells were positive for CD29, CD44, CD73, CD90, CD105, CD133, CD166, HLA class I, Stro-1, and c-Met; 
however, they were negative for CD14, CD31, CD34, CD119, HLA-DR, and c-Kit. Thus, hcMSCs expressed the cell surface antigens typical 
of MSCs. (B) To determine the multilineage differentiation potential of the hcMSCs, the cells were cultured in adipogenic, chondrogenic, 
hepatogenic, and osteogenic medium to induce differentiation. Cytochemical staining showed that the cells differentiated into adipocytes (evaluated 
with red lipid droplets by oil red O), chondrocytes (evaluated with red glycosaminoglycans by safranin O), hepatocytes (evaluated with red 
glycogen deposits by periodic acid-Schiff (PAS)), and osteocytes (evaluated with black calcified nodules by von Kossa), respectively. (C) To examine 
the neural differentiation potential of hcMSCs, the cells were induced for 14 days, and neural marker expression was analyzed by IF staining. 
The differentiated cells exhibited neuronal and glial morphology and expressed several neural markers (GFAP, Tujl, and MAP2), indicating that 
the hcMSC line has neural differentiation potential, as well as adipogenic, chondrogenic, hepatogenic, and osteogenic differentiation potential. 
(D) To confirm the neural differentiation of the hcMSCs at the molecular level, mRNA expression was analyzed. RT-PCR showed that the expression 
of neuronal and glial marker genes, including NF-M, nestin, MAP2, GFAP, and Hes/, increased during neural differentiation. (E) To assess the 
immunosuppressive activity of hcMSCs, lymphocyte proliferation in a mixed lymphocyte reaction was measured by [^H]-thymidine incorporation 
in the presence of hcMSCs. hcMSCs significantly inhibited the proliferation of activated lymphocytes when co-cultured at a ratio of 1:5 
(hcMSCs:PBMCs). (-) Ctrl, negative control; *p = 0.01. 
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were positive for c-MET, but negative for GDI 19 and c-Kit. 
The hcMSCs significantly inhibited T cell proliferation in 
mixed lymphocyte reactions, suggesting that they have in- 



trinsic immunomodulatory activity (Fig. IE). As expected, the 
cells successfully transdifferentiated into adipocytes, chon- 
drocytes, hepatocytes, and osteocytes when induced in vitro 



Table I. Top 30 upreguited or downregulated genes in neural differentalton 
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Table I. Continued 



No. Upregulated genes 


Gene 
symbol 


Genbank 
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Fold- 
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Down regulated genes 
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with adipogenic, chondrogenic, hepatogenic, and osteogenic 
medium, respectively (Fig. IB). Additionally, the cells ex- 
hibited a neural phenotype when culaired in the appropriate 
induction medium. After 14 days of differentiation, the cells 
were fixed and stained for IF analysis. When examined by 
confocal microscopy, the induced cells were positive for the 
neuronal and glial markers Tujl, MAP2, and GFAP (Fig. IC). 
GFAP, neurofilament M (NFM), nestin, MAP2, and Hesl 
mRNA levels increased during neural differentiation (Fig. ID). 
The mRNA levels of these neuroglial maker genes were in- 
creased or maintained for up to 21 days. These results 
showed that the hcMSC line could acquire a neural pheno- 
type after appropriate induction. 

Although neuronal expression profiles of BM-derived MSCs 
have been reported (22,24), to our knowledge, no study has 
assessed a homogeneous MSG population. We chose homo- 
geneous hcMSCs rather than heterogeneous MSG pools in or- 
der to improve the reliability of the gene expression data and 
to eliminate the mixed information obtained from a heteroge- 
neous MSG population. Two independent RNA samples from 
the hcMSG line were isolated before and after neural in- 
duction for 14 days and then subjected to microan'ay analysis. 
The average value of the two independent results was used 
for data analysis. 

The microarray probe set included 27,958 genes. A 2-fold 
increase over expression in control cells was selected as the 
threshold for upregulated genes, and a 2-fold decrease was 
selected as the threshold for downregulated genes. With 
these criteria, we identified 1025 upregulated genes and 1415 
downregulated genes. Table I shows the top 30 upregulated 
genes and the top 30 downregulated genes in differentiated 
MSGs. During neural differentiation, the expression of PCSKl, 
MMPIO, and CYCZi increased 42.3-fold, 30.1-fold, and 24.4- 
fold, respectively. Gonversely, ITGA8, DUSP18, and DKKl 
decreased 23.2-fold, 19.4-fold, and 17.8-fold, respectively. 
The genes were sorted and classified into several groups ac- 
cording to gene ontology. We examined the data with respect 
to neuronal differentiation or nervous system differentiation 
(Supplemental Table). Forty-three genes, including Ereg, neu- 
rotrophic tyrosine kinase receptor 3 (NTrk3), Musashi homo- 
log 1 (MSIl), Sema4G, NDRG2, Aniigol, EphrinBl, Notch3, 
and Pax6 were upregulated, whereas only 3 genes {MyD88, 
GDAP2, and RQCDl) were downregulated. NTrk3, NotchS, 
NDRG2, MSIl, Aniigol, Senia4G, and EphrinBl play im- 
portant roles in neurogenesis. The fact that these upregulated 
genes are involved in neuronal differentiation and nervous 



system development provides evidence of the genuine neuro- 
genic potential of MSGs, refuting some controversial studies 
that suggest a false-positive gain of neural phenotype in MSCs 
treated with neural induction reagents. 

To verify the gene expression profile generated by the mi- 
croarray analysis, genes were randomly selected for evalua- 
tion by RT-PGR (Fig. 2). The expression of PCSK, Ereg, 
MMPIO, Cdxl, WntJa, Sox4, and MEPE increased, whereas 
the expression of Sniad3, SmuifZ, NEDD4L, STATl, DAAM2, 
ITGA6, and DKKl decreased. These results were consistent 
with the microarray data, indicating that our microarray ex- 
periments were highly reliable. 

To generate additional gene expression data, we examined 
epigenetic changes by analyzing promoter methylation and 
compared the results with the gene expression profile. DNA 
methylation is a key epigenetic regulatory mechanism for de- 
velopment and gene expression (25,26) that typically occurs 
at CpG dinucleotides. Clusters of GpG dinucleotides, called 
GpG islands, are usually found in mammalian gene promoters 
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Figure 2. RT-PGR confirmation of the gene expression changes 
observed in the microarray analysis of neurally differentiated hcMSCs. 
After the hcMSCs were differentiated for 14 days, mRNA was 
extracted and subjected to microarray analysis. To confirm the 
microarray data, the expression of randomly selected genes that 
changed 2-fold or more in the microarray analysis was assessed by 
RT-PCR. + denotes increase and — represents decrease. 
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(27). Methylation of the CpG islands in a promoter region 
inhibits gene expression, whereas lack of methylation permits 
gene expression (28). Therefore, examining the CpG methyl- 



ation status of promoters may provide further insight into 
gene expression during the neural differentiation of MSCs. 
We isolated genomic DNA from hcMSCs before and after neu- 



Table II. Promoter methylation patterns of some differentially expressed genes 
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Downregulated & hypermethylated genes 


Gene symbol 


GenBank accession no. 
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Ribonucleotide reductase M2 polypeptide 


RRM2 
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3 
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4 
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Natriuretic peptide receptor C/guanylate cyclase C 
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CDC42 effector protein (Rho GTPase binding) 3 
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7 
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Adenylate kinase 5 
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9 
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Aryl-hydrocarbon receptor nuclear translocator 2 
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12 
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13 


RAB35, member RAS oncogene family 


RAB35 


NM_006861 


14 


Chromosome 14 open reading frame 106 
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ral induction for 14 days and used CpG methylation micro- 
array analysis to determine the genome-wide methylation 
profile. 

In differentiated hcMSCs, 1141 promoters were hyper- 
methylated (5: 2-fold increase) and 1276 promoters were hy- 
pomethylated (^2-fold decrease) when compared with pro- 
moters in control hcMSCs. Among the hypomethylated genes, 
we looked for genes that were upregulated > 2-fold in the 
microarray data. Similarly, among the hypermethylated genes, 
we looked for genes that were downregulated ^2-fold in the 
microarray data. We identified 50 genes whose expression 
was upregulated and whose promoters were hypomethylated, 
along with 68 genes with downregulated expression and hy- 
permethylated promoters (Table II). The genes that were up- 
regulated and hypomethylated included NTrk3, Sox4, SYTIII, 
PABP5, GRIA4, NXPH4, FoxCl, BARHL2, PAPLN, and Pax6. 
Among these genes, NTrk3 and Sox4 belonged to the group 
of 30 genes with the highest expression in differentiated 
hcMSCs, as listed in Table I. Most of genes in Table II are 
involved in neurogenesis or nervous system development, 
confirming the high fidelity of our microarray data. Genes 
that were downregulated and hypermethylated included 
NRXN3, Cdc42EP3, NEDD4L, SIRPa , and DEK. 

DISCUSSION 

Microarray technology enables easy comparison of the differ- 
ential expression of a large number of genes. However, 
false-positive results present an obstacle to microarray data 
analysis. Given that microarray experiments are likely to 
show variation in their results, sample quality is another crit- 
ical factor in the analysis. Some studies have determined the 
gene expression profiles of neural differentiation in MSCs us- 
ing microarray analysis (16,24). The most common method 
for isolating MSCs is FicoU-mediated density gradient cen- 
trifugation. The resulting MSC pool is heterogeneous because 
no additional purification process is used in this method. 
Even the defined sorting method, which uses specific anti- 
bodies, might not completely eliminate heterogeneity in the 
isolated MSC population. We previously developed an easy 
and simple isolation method for obtaining hcMSCs called SCM 
(21). Because no reports have described microarray analysis 
of a homogeneous MSC population, we examined the gene 
expression profile of neurally differentiated hcMSCs. The mi- 
croarray results were confirmed with RT-PCR analysis. In ad- 
dition, we combined the microarray data with independent 



promoter methylation chip array data to generate more accu- 
rate information. 

According to the gene ontology analysis, the expression of 
neuronal-specific genes involved in neuronal differentiation 
and nervous system differentiation increased in neurally dif- 
ferentiated hcMSCs (Supplemental Table). Genes involved in 
synaptic transmission, such as EGR3 (listed in Supplemental 
Table), SYT3 (listed in Table II), nicotinic cholinergic receptor 
5 (not listed, 4.5-fold increase), glycine receptor alpha 3 (not 
listed, 4.4-fold increase), protocadherin beta 13 (not listed, 
3.1-fold increase), neurotensin receptor 1 (not listed, 2,1-fold 
increase), and a voltage-dependent calcium channel (not list- 
ed, 2.1-fold increase), were upregulated. Genes encoding 
channel proteins, including a voltage-gated potassium chan- 
nel, a glycine receptor, glutamate receptors, nicotinic chol- 
inergic receptors, and a voltage-dependent calcium channel, 
were also upregulated (not listed). 

The canonical Wnt signaling pathway plays an important 
role in neurogenesis (29,30). According to our microarray 
analysis, the Wnt signaling molecule Wnt7a was highly upre- 
gulated at the mRNA level, whereas DKKl , a Wnt signaling 
inhibitor, was significantly downregulated 17.8-fold (Table I). 
This is consistent with reports describing the functional role 
of Wnt signaling in neurogenesis (31,32). With regard to Wnt 
signaling molecules, the genes encoding Axin2 and FZDIO 
were upregulated, and those encoding TLE4 and DACTl were 
downregulated in our microarray analysis (not listed). These 
results are consistent with those of a recent study demonstrat- 
ing that FZDIO promotes the development of sensory neurons 
(33). 

One of the interesting gene groups identified in this study 
comprised E3 ubiquitin protein ligases, which catalyze protein 
ubiquitination (34,35). Protein ubiquitination is catalyzed by 
three groups of enzymes, ubiquitin-activating El enzymes, 
ubiquitin-conjugating E2 enzymes, and ubiquitin-ligating E3 
ligases. Functionally, the E3 ligase determines which proteins 
are ubiquitinated. There are two types of E3 ligases: the RING 
type and the HECT type (36,37). In humans, most E3 ligases 
are RING-type enzymes; only 28 ligases are HECT-type en- 
zymes, such as those in the HERC family and NEDD4 family 
(36,37). In neurally differentiated hcMSCs, the gene encoding 
NEDD4L was downregulated and its promoter was hyper- 
methylated (Table II). Other E3 ubiquitin protein ligase genes, 
including SmurfZ (not listed, 4.2-fold decrease), WWP2 (not 
listed, 2.3-fold decrease), and NEDD4 (not listed, 6.4-fold de- 
crease), were also downregulated when assessed by micro- 
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array analysis. Notably, these genes belong to the NEDD4 
family of HECT-type E3 ligases (37). Ubiquitination of RaplB 
by Smurf2 is required for neuronal polarity (38). Downregula- 
tion of Smurf2 expression was confirmed by RT-PCR analysis 
(Fig. 2). Smurfl regulates neurite extension through Rho 
GTPase (39). On the other hand, the expression of other E3 
ligase genes increased, including TRIM50 and RNF25 (data 
not shown). Interestingly, TRIM50 and RNF25 are members 
of the RING-type E3 ligase family. Future studies will de- 
termine whether downregulation of specific HECT-type E3 
ubiquitin ligases, such as NEDD4L and Smurf2, is functionally 
correlated with the acquired neuronal phenotype of MSCs. 
Smurf2 might be a good candidate for future study because 
its downregulated expression was verified by RT-PCR (Fig. 
2). 

Finally, it is notable that TNFAIP6 (l6. 8-fold increase) and 
/ZiiJA/" (16. 5-fold increase) were highly upregulated in neu- 
rally differentiated hcMSCs (Table I). TNFAIP6, also called 
TSG-6, plays a critical role in MSC-mediated immunosup- 
pression (40). TNFAIP6 reduces various inflammatory con- 
ditions, including myocardial infarction, zymosan-induced 
peritonitis, and corneal injury (41-43). ILIRN encodes inter- 
leukin-1 receptor antagonist (IL-lRa), which can modulate 
various IL-l-related inflammatory responses by inhibiting IL-1 
a and IL-1 /3 (44). Anakinra, an approved IL-lRa drug used 
to treat rheumatoid arthritis, has therapeutic effects on a varie- 
ty of local, systemic, hereditary autoimmune inflammatory 
diseases (45). 

We evaluated the gene expression profile of neurally differ- 
entiated MSCs using homogeneous clonal cells. To verify the 
gene expression profile, we confirmed the results from the 
microarray analysis with RT-PCR and analyzed data obtained 
from CpG methylation microarray analysis. Our results pro- 
vide evidence for the neuronal and glial potency of MSCs. 
Thus, MSCs can be used for the development of cell therapies 
that target neuroimmune disorders such as stroke, multiple 
sclerosis, and Parkinson's disease. The upregulation of some 
anti-inflammatory genes, such as TNFAIP6 and KIRN, during 
neural differentiation improves the outlook on the therapeutic 
utility of MSCs. Furthermore, neurogenically induced or neu- 
rogenically primed hcMSCs retained their immunosuppressive 
activity against activated T cells in vitro (unpublished data). 
Thus, an optimized balance between the immunomodulatory 
activity and neurogenic potency of MSCs may provide a new 
therapeutic solution for incurable neuroimmune diseases. 
Although our study is only one small step toward under- 



standing the biological and molecular properties of MSCs, fur- 
ther investigation should contribute to the development of 
stem cell therapy. 
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